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Abstract. The most detailed calculations of the p- 
process call for its development in the /Ne layers of Type 
II supernovae. In spite of their overall success in reproduc- 
ing the solar system content of p-nuclides, they suggest a 
significant underproduction of the light Mo and Ru iso- 
topes. On grounds of a model for the explosion of a 25 
Mq star with solar metallicity, we demonstrate that this 
failure might just be related to the uncertainties left in the 
rate of the 22 Ne (a , n) 25 Mg neutron producing reaction. 
The latter indeed have a direct impact on the distribution 
of the s-process seeds for the p-process. 



1. Introduction 

The most successful p-process models available to-date 
call for the synthesis of the stable neutron-deficient nu- 
clides heavier than Fe in the O/Ne layers of Type II su- 
pernovae (SNII) (Rayet et al. 1995, hereafter RAHPN). 
In spite of their many virtues in reproducing the solar- 
system p-nuclide abundance distribution, they however 
suffer from some shortcomings. One of them concerns their 
persistent underproduction of the light Mo ( 92 Mo, 94 Mo) 
and Ru ( 96 Ru, 98 Ru) p-isotopes. Some have tried to rem- 
edy this situation with exotic solutions, calling in partic- 
ular for accreting neutron stars or black holes (e.g. Schatz 
et al. 1998). The level of the contribution of such sites to 
the solar system content of the nuclides of concern here is 
quite impossible to assess in any reliable way. In contrast, 
it has been emphasized many times over the last decade 
that the problem might just be due to some misrepresen- 
tation of the production in the He-burning core of massive 
stars of the s-nuclide seeds for the p-process (e.g. Arnould 
et al. 1998). 

The aim of this Letter is to scrutinize the latter, 'non- 
exotic', solution in a quantitative way by duly taking into 
account the uncertainties that still affect the rate of the 
22 Ne (a , n) 25 Mg reaction, as they appear in the NACRE 
compilation of reaction rates (Angulo et al. 1999). Clearly, 
these uncertainties in the key neutron producer in condi- 
tions obtained during central He burning in massive stars 



have a direct impact on the predicted abundances of the s- 
nuclide seeds for the p-process, as already analyzed quan- 
titatively by Meynet & Arnould (1993). Another poten- 
tial embarassment of the p-process predictions identified 
by RAHPN is a SNII overproduction of oxygen relative to 
the p-nuclides. We show that this problem might be cured 
along with the one of the underproduction of the light Mo 
and Ru isotopes if the 22 Ne (a , n) 25 Mg rate is modified 
adequately within a range permitted by the NACRE com- 
pilation. For the sake of illustration, we just discuss here 
the case of a 25 M© solar metallicity (Z = Zq) star. A 
more complete study dealing in particular with a set of 
stars with different masses and metallicities, and analyz- 
ing the impact of the uncertainties in the rates of a variety 
of reactions, is currently under way. 

The adopted input physics is briefly described in Sect. 
2, and some results are presented in Sect. 3. Conclusions 
are drawn in Sect. 4. 



2. Input physics 

The calculations reported here are based on a model for 
a Zq 8 Mq helium star already considered by RAHPN. 
It corresponds to a main sequence mass of about 25 M© 
and is evolved from the beginning of core helium burn- 
ing to the supernova explosion. Details about this model 
can be found in Hashimoto (1995), and are summarized in 
RAHPN. As in RAHPN, 20 O/Ne-rich layers with explo- 
sion temperatures peaking in the (1.8-3.3) x 10 9 K range 
are selected as the P-Process Layers (PPLs). Their total 
mass is approximately 0.58 M©. The deepest PPL is lo- 
cated at a mass of about 1.94 M©, which is far enough 
from the mass cut for all the nuclides produced in this 
region to be ejected during the explosion. 

The p-process reaction network and its numerical 
solver are described by RAHPN. A series of their selected 
nuclear reaction rates are updated, however. In particular, 
the NACRE 'adopted' rates are used for charged particle 
captures by nuclei up to 28 Si. For heavier targets, the rates 
predicted by the Hauser-Feshbach code MOST (Goriely 
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1997) are used, except for the experimentally-based neu- 
tron capture rates provided by Beer et al. (1992).^ 

As already pointed out in Sect. 1, we turn our special 
attention to the impact of the uncertainties remaining in 
the rate of 22 Ne (a , n) 25 Mg . For temperatures of about 
2 — 3 x 10 8 K at which the s-process typically develops 
during core He burning in massive stars (e.g. Rayet & 
Hashimoto 2000), the NACRE upper limit on this rate 
is 50-500 times larger than the 'adopted' value (see An- 
gulo et al. 1999 for details). In order to quantify the con- 
sequences of this situation for the predicted abundance 
distribution of the s-nuclide seeds for the p-process, and 
ultimately for the p-nuclide yields themselves, we perform 
nucleosynthesis calculations for five different rates rang- 
ing from the NACRE adopted value to its upper limit. 
These rates, labelled Ri (i=l to 5) in the following, are 
defined and displayed in Fig. 1. Q They are used in the Z Q 
25 M Q star referred to above to calculate the abundances 
of the s-process nuclides at the end of core He burning. 
The results are shown in Fig. 2 for the s-only nuclides. 
Use of Ri leads to the classical 'weak' s-process compo- 
nent pattern (e.g. Rayet & Hashimoto 2000), exhibiting 
a decrease of the overproduction (with respect to solar) 
of the s-nuclides by a factor ranging from wlOO to about 
unity when the mass number A increases from about 70 to 
100. In the heavier mass range, the s-process 'main com- 
ponent' supposed to originate from low- or intermediate- 
mass stars takes over. This 'canonical' picture changes 
gradually with an increase of the 22 Ne (a , n) 25 Mg rate, 
more 22 Ne having time to burn, releasing more neutrons, 
before He exhaustion in the core. The direct result of this 
is a steady increase of the overproduction of heavier and 
heavier s-nuclides. For example, with the extreme R5 rate, 
the overproduction factor increases from 10 3 to 10 4 for A 
varying from about 70 to 90, before decreasing to a value 
around unity for A « 150 only. 

At first sight, it might be felt that the s-process abun- 
dance distributions obtained with large enough Rj values 
exhibit some unwanted or embarrassing features. One of 
these concerns the underproduction of the A w 70 — 76 s- 
nuclides relative to the A « 80 — 90 ones. Another one re- 
lates to the fact that a more or less substantial production 
of heavy s-nuclides (like in the Ba region) would screw up 
the pattern of the s-process main component ascribed to 
lower-mass stars. In our opinion, none of these predictions 
can really act as a deterrent to 22 Ne (a , n) 25 Mg rates sub- 
stantially in excess of Ri . On the one hand, the absence of 
ab initio self-consistent calculations of the s-process in low- 



The NACRE and MOST rates are available in the Brussels 



Nuclear Astrophysics Library (http://www-astro.ulb.ac.b6) 

2 The calculations reported here were completed when we 
have been informed of new low-energy measurements of the 
22 Ne (a , n) 25 Mg cross section (J.W. Hammer & M. Jaeger, 
private communication). Rate estimates based on these new 
data are not available yet. It seems, however, that the revised 
upper limit might range somewhere between R4 and R 5 
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Fig. 1. The five 22 Ne (a , n) 25 Mg rates Ri used in our 
calculations for temperatures of relevance for core He 
burning in massive stars. All rates are normalized to Ri, 
which is the NACRE 'adopted' rate. R5 corresponds to 
the NACRE upper values. R3, R2 and R4 are the geomet- 
rical means between Ri and R5, Ri and R3, and R3 and 
R5, respectively (see also footnote 2) 
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Fig. 2. Distribution of the abundances, normalized to so- 
lar values, of the s-only nuclides at the end of core He 
burning in the considered Zq 25 M© model star, for the 
22 Ne (a , n) 25 Mg rates R, (i = 1 to 5) denned in Fig. 1, all 
the other ingredients of the model being kept unchanged 
(see Rayet & Hashimoto 2000). The dashed line is the 
distribution adopted by RAHPN for their p-process cal- 
culations 



and intermediate mass stars does not allow at this time to 
predict the exact shape of the main s-process component 
which is classically assigned to these stars. As a conse- 
quence, a contribution to the main component by massive 
stars cannot be excluded, even if it may disturb some tra- 
ditional views on the subject. On the other hand, the re- 
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duction of the light s-process nuclide production by mas- 
sive stars could well be compensated by their increased 
synthesis by some low- or intermediate-mass stars when 
rates larger than Ri are considered (Goriely & Mowlavi 
2000). The classical 80 Kr overproduction problem found in 
the massive star s-process (e.g. Rayet & Hashimoto 2000) 
could also be eased with increased 22 Ne (a , n) 25 Mg rates, 
as demonstrated by Fig. 2. For these same rates, note that 
80 Kr is not overproduced either in some of the calculations 
of Goriely & Mowlavi (2000) which predict high yields of 
the other light s-nuclides. 

Figure 2 also suggests that a discrepancy, if any, be- 
tween the observed Ba overabundance in the SN1987A 
ejecta and the model predictions could be cured in a nat- 
ural way by increasing the adopted 22 Ne (a , n) 25 Mg rate. 
The [Ba/Fe]sAr/[Ba/Fc]£MC ratio is observationally still 
quite uncertain, values between about 5 and 20 having 
been reported (e.g. Mazzali & Chugai 1995). Prantzos et 
al. (1988) have calculated lower values of 2.6 to 4.7 with 
the 22 Nc (a , n) 25 Mg rate of Fowler et al. (1975). This rate 
is on average comparable to Ri in the temperature range 
of relevance to the s-process. We have not conducted any 
new s-process calculation in a specific SN1987A progenitor 
model. Instead, some rough estimates based on the pro- 
cedure of Prantzos et al. (1988) in which their adopted s- 
process Ba mass fraction is replaced by the one calculated 
for the model star adopted here have been made. Assum- 
ing that the LMC metallicity is one third of the solar one, 
we predict [Ba/Fe]sAr/[Ba/Fe]LMC ratios from 3 to 14 for 
rates increasing from Ri to R5. Theory could thus ac- 
count for quite substantial SN1987A Ba productions with 
high enough 22 Ne (a , n) 25 Mg rates (compatible with the 
NACRE data). 

As discussed by RAHPN, it is a fair approximation 
to adopt the s-process abundance distributions of Fig. 2 
as seeds for the p-process. For the A < 40 species, the 
initial abundances in the PPLs are taken from the detailed 
stellar models. Although these models have been obtained 
with rates that may differ from the NACRE ones adopted 
here, this inconsistency is certainly not responsible for any 
intolerable distorsion in the predicted s-process seeds or p- 
process yields. 

3. Results and discussion 

The various seed abundances of Fig. 2 are used to compute 
the production of the p-nuclidcs in the PPLs of the Z Q 25 
Mq star considered here. As in RAHPN, the abundance of 
a p-nuclide i is characterized by its mean overproduction 
factor (Fi) = (JQ) /A iiQ , where Xj iQ is its solar mass 
fraction (Anders & Grevesse 1989), and 
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J2( x i,n + X iin -i)(M n - M„_i)/2, 
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where X i n is the mass fraction of isotope i at the mass 
coordinate M n , M p = J2 n>1 (M n — M n _i) is the total 
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Fig. 3. Values of (Fi) /F and of F derived from the seed 
abundances calculated with the 22 Ne (a , n) 25 Mg rate Ri 
(upper panel), R3 (middle panel) and R5 (lower panel). 
Lines connect different isotopes of the same element. The 
dotted horizontal lines delineate the 0.3 < (Fi) / F < 3 
range 

mass of the PPLs, the sum running over all the PPLs 
(Mo corresponds to the bottom layer). An overproduc- 
tion factor averaged over all 35 p-nuclei is calculated as 
Fq = (Fi) /35, and is a measure of the global p-nuclide 
enrichment in the PPLs. So, if the computed p-nuclei 
abundance distribution were exactly solar, the normalized 
mean overproduction factor (Fi) / F would be equal to 
unity for all i. 

Figure 3 shows the normalised p-nuclide overproduc- 
tion factors derived from the seed abundance distributions 
calculated with the 22 Ne (a , n) 25 Mg rates Ri, R 3 and R 5 . 
Changes in the shape of the p-nuclide abundance distribu- 
tion are clearly noticeable, at least for A ^ 120. The use of 
Ri leads to a more or less substantial underproduction of 
not only 92 Mo, 94 Mo, 96 Ru and 98 Ru, a 'classical' result in 
p-process studies (see RAHPN), but also of 78 Kr and 84 Sr, 
which was not predicted in previous calculations. This new 
feature directly relates to the larger abundances around 
A « 80 used by RAHPN (dashed curve in Fig. 2), in con- 
trast to the much flatter seed distribution obtained with 
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Fig. 4. Values of the normalized overproduction factors 
of the Kr, Sr, Mo and Ru p-isotopes as a function of the 
22 Ne (a , n) 25 Mg rate 

Ri. This Kr-Sr-Mo-Ru trough is gradually reduced, and 
in fact essentially disappears, for 22 Ne (a , n) 25 Mg rates 
of the order or in excess of R3. This situation is most 
clearly illustrated by Fig. 4. In these very same conditions, 
(Fi) /Fq for 113 In and 115 Sn comes much closer to unity as 
well. It has to be noticed that this situation does not result 
from a stronger production of these two nuclides by the 
p-process, but instead from their increased initial abun- 
dances associated with a more efficient s-process when 
going from Ri to R 5 . In contrast, the (Fi) /F pattern 
does not depend on the adopted 22 Ne (a , n) 25 Mg rate for 
A <; 140. This is expected from a mere inspection of the 
s-nuclide seed distributions displayed in Fig. 2. In partic- 
ular, 152 Gd and 164 Er remain underproduced. This cannot 
be considered as an embarassment as these two nuclides 
can emerge from the s-process in low- or intermediate- 
mass stars. 

In addition, the overall efficiency of the p-nuclide pro- 
duction substantially increases with increasing 22 Ne burn- 
ing rates. More specifically, F is multiplied by a factor of 
about 15 when going from Ri to R 5 . This could largely 
ease, and even solve, the problem of the relative underpro- 
duction of the p-nuclides with respect to oxygen identified 
by RAHPN. For their considered 25 M Q model star cal- 
culated with the 12 C (a , 7) 16 rate from Caughlan et al. 
(1985), they obtain F = 130 and report a value of 4.4 
for the ratio of the oxygen to p-process yields. This value 
would come close to unity for 22 Nc (a , n) 25 Mg rates in 
the vicinity of R3-R4, as the p-nuclides would be about 3 
to 6 times more produced than in RAHPN. 

4. Conclusions 

This Letter makes plausible that the long-standing puz- 
zle of the underproduction with respect to solar of the 



p-isotopes of Mo and Ru in SNII explosions could be 
quite naturally solved by just assuming an increase of the 
22 Ne (a , n) 25 Mg rate over its 'nominal' value. More specif- 
ically, this could be achieved by multiplying the NACRE 
'adopted' rate by factors of about 10 to 50 in the tem- 
perature range at which the s-process typically develops 
during core He burning in massive stars. These factors are 
well within the uncertainties reported by NACRE. As an 
important bonus, this increased rate would also largely 
avoid (i) the underproduction of 78 Kr and of 84 Sr which 
we predict here for the first time to be concomitant to 
the light Mo and Ru one, (ii) the too low production of 
113 In and 1 Sn, and (iii) the overall underproduction of 
the p-nuclides with respect to oxygen noted by RAHPN. 
In direct relation with an increased 22 Ne (a , n) 25 Mg rate, 
more s-process Ba could also be ejected by SNII events. 
Our predictions confortably overlap the range of Ba over- 
abundances reported for SN1987. 

This array of pleasing features has of course not to be 
viewed as a proof of the validity of the assumption that the 
true 22 Ne (a , n) 25 Mg rate is higher than usually thought. 
It may just be a hint that there might be ways around ex- 
otic solutions. This conclusion applies at least if one relies 
on the simplistic (and the only ones to be available for 
our purpose) supernova models used here and in previous 
p-process calculations (see RAHPN et references therein), 
as well as in a myriad of other explosive nucleosynthesis 
calculations. 
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